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I review the state of the art of the investigation on the structure formation in /(_R)-gravity based 
on the Covariant and Gauge Invariant approach to perturbations. A critical analysis of the results, 
in particular the presence of characteristic signature of these models, together with their meaning 
and their implication is given. 



' I. INTRODUCTION 

£SJ , Determining the nature of dark energy [l[ has become, undoubtedly, one of the central issues in modern cosmology 
^ ■ as well as in theoretical physics. The interest in this problem has grown to such level that nowadays much of the 
\ research work in theoretical cosmology is focused on the quest of a theoretical framework in which dark energy can 
["T | ■ be explained. One idea that has recently gained much popularity is that dark energy has a geometrical nature i.e. it 
' is originated not by some unknown new form of energy density that dominate the cosmic evolution, but by the fact 
. that the gravitational interaction, at least on cosmological scales has to be modified. 

Among the many modified versions of General Relativity (GR), higher order theories of gravity, and particularly 
r— 1| one of its subclass, called /(i?)-gravity, has been thoroughly studied (see e.g. for some reviews). In these models 

^ i the action for the gravitational interaction is written as a generic analytic function of the Ricci scalar and they can 
I | be recovered naturally in the low energy limit of fundamental theories (see e.g. [f|). However, the reason why those 
H„> models have attracted such attention is that it has been proven in many different ways that these theories are able to 
l_j reproduce in a natural way the cosmological "footprint" of Dark Energy i.e. the cosmic acceleration [1, [H, @-[8[ . Such 
phase cannot be achieved if one considers only GR without the introduction of some additional field with suitable 
(and, unfortunately, odd) thermodynamical properties. 

However, /(i?)-gravity is not the only model capable to achieve this kind of results and the modification of a theory 
as well established as GR requires a serious, careful and unbiased analysis. To be completely fair so far although 
£q \ many concerns have been raised, there is no definite argument against /(i?)-gravity or, for what matters, an argument 
■ that proves definitively the presence of corrections to the gravitational interaction. Therefore there is a great need to 
devise observable features of these models and compare them to the observations. In this respect cosmology is of great 
help. In fact, the cosmological models of any theory of the gravitational interaction are by far the easiest application 
' of these schemes and, because of the very nature of the cosmological processes, their analysis is able to give insights 
—i , of the behavior of these theories in a wide range of energy scales. 

Following this spirit in this paper I will review what it has been discovered so far on the behavior of the scalar 
. ! (matter) perturbations and the structure formation in this framework. On this topic many papers have been published 
that employ a variety of different techniques/approximations In what follows I will only consider the results 
obtained with a specific method: the Covariant Gauge Invariant (CoGI) approach. This approach was proposed 
at the beginning of the eighties by Bruni and Ellis [l(J and generalized also by Hwang and Dunsby TJl4|. It is 
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constructed on a formalism optimized to analyze cosmological models: the 1+3 covariant approach THT" 16|. Using 
the CoGI approach we will realize that the behavior of the perturbation in this type of fourth order gravity can be 
very different from the one of the standard picture, but not necessarily incompatible with the observations. I will 
show that the process of structure formation offers in some cases very specific signatures, which can be easily tested 
with the data currently available. 

The paper will be organized as follows. In Section II the general equations for /(i?)-gravity are briefly described. 
In Section III I will review briefly the 1+3 covariant approach to cosmology and its application to /(-R)-gravity. In 
Section IV on the base of the 1+3 approach the CoGI approach is developed for the study of the scalar perturbations. 
In Section V I will summarize the results obtained for some simple specific models. Finally section VI is dedicated to 
the conclusions. 

Unless otherwise specified, natural units (K = c = ks = SttG = 1) will be used throughout the paper, Latin 
indices running from to 3. The symbol V represents the usual covariant derivative and d corresponds to partial 
differentiation. I use the — , +, +, + signature and the Riemann tensor is defined by 

R a bcd = W a bdiC - W a bc , d + W e bd W a ce - W f bc W a df , (1) 
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where the W a bd are the Christoffcl symbols (symmetric in the lower indices), defined by 

1 



in- 



fed 



:9 ae (9be,d + 9ed,b - gbd,e) ■ 



The Ricci tensor is obtained by contracting the first and the third indices 

Rab = g Cd Racbd ■ 

Symmctrization and antisymmctrization over the indices of a tensor are defined as 



T(ab) = ^ (Tab + Tba) 



T[ab\ = ^ (Tab - Tba) 



Finally, the Hilbert-Einstein action in the presence of matter is given by 

A= J dx 4 y/=g[R + 2£ m ] . 



(2) 
(3) 

(4) 
(5) 



II. GENERAL EQUATIONS FOR FOURTH ORDER GRAVITY. 



In four dimensional homogeneous and isotropic spacetimes i.e. Friedmann Lemaitre Robertson Walker (FLRW) 
universes, a general action for fourth order gravity can be written as an analytic function of the Ricci scalar only: 



A = 



-g [f{R] + 2jC Ti 



(0) 



where C m represents the matter contribution. Varying the action with respect to the metric gives the generalization 
of the Einstein equations: 



Rab - \ 9abR) = T2 + l -g ah (/ - Rf) + V b V a /' - ffa6 V c V c /' , 



(7) 



where / = f(R), f 



df(R) 



and T, 



2 6(y/=gC m ) 



represents the stress energy tensor of standard matter. 



dr y/=g 5g a b 

These equations reduce to the standard Einstein field equations when f(R) = R. It is crucial for our purposes to be 
able to write ((7]) in the form 1 



rpVa i rpli rptot 

ab ab < ab ab ' 



(8) 



where IT 



T, 



r 



and 



rpR _ J_ 

ab fl 



^Jab (/ - Rf) + VfeV a /' - <?ahV c V C /' 



(9) 



represent two effective "fluids" : the curvature "fluid" (associated with X^) and the effective matter "fluid" (associated 
with f%l) [H HH. In this way fourth order gravity can be treated as standard Einstein gravity plus two "effective" 
fluids and we can adapt easily many of the techniques developed for GR-based models. 

Let us look at the conservation properties of these effective fluids using the Bianchi identities Q ■ The covariant 
derivative of the total stress energy momentum in (j8]) yields: 



= \7 b T^ = 



ab 



f ^T™ b V b R + V b T* 



(/') 



(10) 



1 For this step to make sense it is crucial that we suppose f'(R) ^ at all time. This could be problematic when one deals with cosmologies 
in which the Ricci scalar is zero (like the radiation dominated ones) and functions / such that /'(0) = 0. In the following, although we 
will retain a general barotropic factor, we will consider matter to be dust in all but one case. 
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Using the field equations and the definition of the Riemann tensor, it is easy to show that the sum of the first two 
terms on the RHS of the previous expression is zero. Thus, T*°^ ' b oc Tj^J . b and the total conservation equation reduces 
to the one for matter only. A general, independent proof of this result has been given by Eddington [TtJ and then by 
others (like in [lj|). They showed that the first variation for the gravitational action is divergence free regardless of 
the form of the invariants that we choose for the Lagrangian. This means that no matter how complicated the effective 
stress-energy tensor T'°£ is, it will always be divergence free if T£ b l . b = 0. As a consequence, the total conservation 
equation reduces to the one for standard matter only. For our purposes this is important because it tells us that no 
matter how the effective fluids behave, standard matter still follows the usual conservation equations T^J ' = 0. 

The form @ of the field equations allows us to use directly the 1+3 covariant approach that will be introduced in 
the next section. 



III. THE 1+3 COVARIANT APPROACH TO COSMOLOGY 



Basically all the calculations and the results that will be given in this paper are based on the 1+3 covariant 
approach [15l. HtSj. Such approach is, conceptually, not different from the ADM one, but it is specifically adapted to 
treat homogeneous and isotropic spacetimes and can simplify quite a lot the calculations. Therefore, before starting, 
it is worth to give a very brief review of this approach and its application to /(i?)-gravity 

The 1+3 approach is based on the choice of a specific family of preferred worldlines which represent specific classes 
of observers and can be associated to a timelike vector field u a . Using this field one can split the metric tensor as 

9ab = h a b - U a Ub , (11) 

i.e. the spacetime is foliated in hypersurfaces with metric h a b orthogonal to the vector field u a . In this way any 
affinc parameter on the worldlines associated to u a can be chosen to represent "time" and the tensor h a b (h a c h c b = 
h a b , h a a — 3 , h a b u b = 0) determines the geometry of the instantaneous rest-spaces of the observers we have chosen. 
Using u a and h a b, one can then define the projected volume form r\ a b c = u d r\ a bcd on the rest spaces, the covariant time 
derivative (') along the fundamental worldlines, and the fully orthogonally projected covariant derivative V: 

X ab cd = u e W e X ab cd , V e X ab cd = h a f h b g h p c h q d h r e V r X f9 pq . (12) 
Also, performing a split of the first covariant derivative of u a into its irreducible parts, namely 

V a u b = - u a a b + - 6 h ab + a ab + u ab , (13) 

one can define, in analogy with classical hydrodynamics [Tsl. [20l|. the basic kinematical quantities of this formalism 
[Til ]. In (|13[) . = V a u a represents the rate of volume expansion of the worldlines of u a so that the standard Hubble 
parameter is H: H = 30; a ab = V ( a u b) is the trace-free symmetric rate of shear tensor describing the rate of 

distortion of the observer flow and we have o a b = 0Y O &); CabU 1 ' = 0, cr a a = 0; uj a b = V[ a u&] is the skew-symmetric 
vorticity tensor describing the rotation of the observers relative to a Fermi-propagated (non-rotating) frame and 
Lu a b = U3[ a b], ^>abU b — 0; ab — Ub is the acceleration vector, which describes the non-gravitational forces acting on the 
observers. 

The matter energy- momentum tensor T a b of a general fluid can also be decomposed locally using u a and h a b'- 

Tab = V u a u b + q a Ub + u a qb+P h ab + n ab , (14) 

where u = (T a bU a u b ) is the relativistic energy density relative to u a , q a = —Tb c u b h ca is the relativistic momentum 
density, which is also the energy flux relative to u a (q a u a = 0), p = \(T a bh ab ) is the isotropic pressure, and 
T^ab = T cd h c ( a h d b) is the trace- free anisotropic pressure for which 7r a a = , ir a b — ""(afc)- This allows us to extract 
information on the thermodynamics associated to this fluid. 

The kinematics and thermodynamics quantities presented above completely determine a cosmological model. The 
advantages in using these variables is that they allow a treatment of cosmology that is both mathematically rigorous 
and physically meaningful and they are particularly useful in the construction of the theory of perturbations. Their 
evolution and constraint equations, also known as 1+3 covariant equations (see [HI), are completely equivalent to the 
Einstein equations and characterize the full evolution of the cosmology. 
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A. The 1+3 covariant approach for /(i?)-gravity 



Following the above scheme let us apply the 1+3 formalism to /(i?)- gravity. As first step one needs to choose 
suitable frame, i.e., a 4-velocity field u a . Following [l_| ___]], we will choose the frame it™ 1 comoving with standard 
matter represented by galaxies and clusters of galaxies. This frame basically coincides to our specific point of view: 
as earth bound observers we are comoving with these object. We will also assume that in u™ standard matter is a 
barotropic perfect fluid with equation of state p m = wfi m . 



Relative to it m , the stress energy tensor T*£* given in © can be decomposed as 

w 1 



/' 



1 ab u u — P T P ■ 



V 



totub„,c _ ~m , R 



T tot h b n 
1 bc n a u 



^h ab =p m + P R 

rptot lc -ud ~ m I 

- J cd n <a a b> — n ab "+ 



j 6 ? 



with 



l ab 



n ab 

r 



(15) 
(16) 

(17) 



Since we assume that standard matter is a perfect fluid in it™, q™ and are zero, so that the last two quantities 
above also vanish. The effective thermodynamical quantities for the curvature "fluid" are 



1 

7 
i 

7 

2 



(Rf - /) - ef'R + f'V 2 R + /"' V 6 i?V 6 i? 



1 

2 

^(/ - Rf) + f'R + f"'R 2 + l&f'R lf"V 2 R+ 



f"V a RV a R + f'a b V"R 



9? 



1 

'7 
i 
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f"RV a R + f'VaR - 2®f"VaR 



f"V {a V b) R + f"V {a RV b) R - a Qb /"i? 



(18) 

(19) 
(20) 
(21) 



The twice-contracted Bianchi Identities lead to evolution equations for /i m , fi R , q R and are given in [T_|. In this way 
the 1+3 equations for /(i?)-gravity in the frame it m can be written as: 



+ \<d 2 + (J ab a ab - 2u a io a - V a u a + u a u a + \{p" 1 + 3p m ) = + 3p fl ) 



Expansion propagation (generalized Raychaudhuri equation): 

i 

Vorticity propagation 
Shear propagation: 



(j ( a ) + |6w a + \cw\u a - a ab uj b = 



°(ab) + §©Cab + E ab - V( a 1t 6 ) + cr c ( Q crt) C + W( Q a; b ) - U( a U b ) = \^ab~ 

Gravito-electric propagation: 

E {ab) + OE ab - cm\H ab + i(/i m + p m )<r ab - 2ii c r, cd(a H b) d - 3<r c{a E b) c + uj c T] cd[a E b) d 

= -llM Mob _ 2 V (ab) - 2 V (a1b) - 6 W7r a6 - 2° (a V b)c ~ 2 U Vc{a^b)d ■ 

Gravito-magnetic propagation: 

H {ab) + GH ab + cwclE ab - 3a c{a H b) c + uj c ri cd(a H b) d + 2ii c ri cd(a E b) d 

R I 1 c ,_, di? 



Vorticity constraint: 



icurl7rf b - %u)( a q'b) + %<r"(aVb)cQd 



V a w a - ii a Lu a = 



(22) 
(23) 
(24) 

(25) 

(26) 
(27) 
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Shear constraint: 

V b a a6 - curiae - §V a 9 + 2[w,u] a = -g? . (28) 

Gravito-magnetic constraint: 

curlcr afc + V( a a; 6 ) - i? afc + 2?i (a u; b ) = . (29) 

Gravito-electric divergence: 

V b E ab - ±V a M m - [a, H] a + 3H ab oj b = \a b a q R - g fl ] Q - |V 6 7rf 6 + §V a A* fl - |0gf . (30) 

Gravito-magnetic divergence: 

V b H ab - (fl m + p m )u a + [a,E] a - 3E ab cj b = -£curlg* + (fi R + p R )w a - \\a^ R \ a - \i^J> . (31) 

Standard Matter Conservation (twice contracted Bianchi identities) 

(T = -Q(n m +p m ), (32) 
V°p m = -(p m +p m )u a . (33) 

Curvature fluid Conservation (twice contracted Bianchi identities) 

lH R + V a 9 f = - 6 (u." +p R )-2 (u a q R ) - (a^ R a ) +H m ^, (34) 
Q(a) + V aP R + V b n R =-lGq R - a b q R - + p R ) u a - u b n R - V bc u b q R + fj, m (35) 

In the equations above the spatial curl of a vector and a tensor is 

(curl X) a = rf bc V b X c , (curl X) ab = r? cd(o VcX h) d , (36) 

respectively. 

Finally, ui a = ^Va bc u bc and the covariant commutators are 

[X, Y] a = i lacd X c Y d , [W, Z] a = Vacd W c e Z de . 

Note that for f(R) = R, one has f'(R) = 1 and fJ, R ,p , q R , ir R b = and therefore the above equations reduce to the 
ones for GR. Using (|22 p -(|35 [l one can analyze any type of f(R) cosmology. Such analysis can be performed either 
directly or with the use of alternative techniques, like the Dynamical System Approach (see e.g. @, H HE HI] ) ■ 

IV. COVARIANT GAUGE INVARIANT PERTURBATIONS 

Once the evolution of a cosmological model has been determined, one can investigate the perturbations around a 
specific metric which is a solution for the background cosmology 2 . In order to analyze the behavior of the perturbations 
we will use the 1+3 covariant approach described above to construct a Covariant Gauge Invariant (CoGI) theory of 
perturbations. Such theory will characterize the evolution of the perturbations as covariant equations written in terms 
of variables which have a very straightforward physical meaning and are naturally gauge invariant. 

The first step is the choice of the background. In the CoGI approach this is not done by assigning a metric, 
but rather by recognizing which 1+3 quantities are zero in the background and which are not. In what follows we 
will consider expanding (O =/= 0) homogeneous and isotropic (a ab = 0, u> ab = 0) backgrounds. In this setting we 
will characterize the perturbations in terms the of 1+3 quantities seen in the previous Section and their projected 
gradients. 



2 It is worth to stress here that this last prescription is a particularly important one. Choosing a random background which is incompatible 
with the order cosmological equations can potentially destroy the entire predictive power of the perturbation equations. 
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For example, to characterize the energy density perturbation, the natural variable is X a = V a /i. This vector 
represents any spatial variation of the energy density fi (i.e. any over-density or void) and it is in principle directly 
measurable pi|. However, a more suitable quantity to describe density perturbation is 

D a = -X a = -V aM , (37) 
M I 1 

where the ratio X a j \x allows one to evaluate the magnitude of density perturbations relative to the background energy 
density and the presence of the scale factor S guarantees that it is dimensionless and comoving in character. The 
magnitude of D a (D = ^/D a D a ) is closely related with the quantity Sfj,/fj,, but represents a real covariant and gauge 
invariant spatial fluctuation. In fact, it is easy to show that the Bardeen variable e m which corresponds to Sfi/fj, in 
the comoving gauge, is the scalar harmonic component of D a (i.e. its scalar "potential" ) (l4| . 
Other important quantities relevant to the evolution of density perturbations in GR are 

Z a ee SV a e , C a ee S 3 V a R 3 , (38) 

which represent, the spatial gradient of the expansion and the spatial gradient of the 3-Ricci scalar respectively. These 
variables are not independent from Q37p , but they are related by a constraint coming from the spatial derivative of the 
Gauss equation 0, [2l|. Moreover it can be proven that these variables, as well as any other quantity which vanish 
in the background, are gauge- invariant (25[. 

The variables defined above contain a lot of information, such as evolution of non-spherically symmetric structures 
and vortical motions, which is not directly associated to the evolution of matter fluctuations. In the following we will 
focus on this last type of perturbations only and, specifically, on the ones directly related to spherically symmetric 
collapse. To extract this information from the variables (|37[) and (|38|) . we use the local splitting 



V a X b — X ab — -h ab X + + X[ ab ] , (39) 

where 

X = V a X a , (40) 



and 



^ab — X (ab) - -^h ab X . (41) 



Because of its geometrical properties it is clear that only the scalar part X can describe a spherically symmetric 
collapse. As a consequence, from now on, we will deal only with the scalar part of the perturbation variables. Such 
part can be extracted applying the comoving differential operator SV a to (p7)l and (|38l) : 



S 2 ~ 

A™ = — VV n , Z = 5 2 V 2 6 , C = S A W 2 R 3 ■ (42) 

Clearly these variables are gauge invariant, for the same reasons for which D a , Z a , C a are. In dealing with /(i?)-gravity, 
however, (|4"2"j) will not be able to describe all the degrees of freedom of the theory. In the next section we will see how 
this can be addressed. 



A. CoGI Perturbations for /(i?)-gravity 

Let us now derive explicitly the perturbation equations around an homogeneous and isotropic background for f(R)- 
gravity as seen by an observer associated to u™ and in presence of a fluid which is perfect in this frame and has an 
equation of state p m = wfj, m . In this setting the (f22|) - (|35l) give the following background equations: 

e 2 = 3^ + 3 M * - , (43a) 

e + ie 2 + ±j( P m + 3 P m ) + ±( P R + ip R ) = o , (43b) 

/i m + 8(Ai m +p m ) = 0, (43c) 

ell 
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where /j, R and p R are given in (fT8|) and (|T9l) and A3 = 6K/S 2 with the spatial curvature index A = 0, ±1. 

As mentioned above /(iJ)-gravity contains more degrees of freedom than standard GR. Therefore in order to 
complete the description of the evolution of the perturbations in this framework we will need to define some additional 
perturbation variables. Following the tradition to consider the Ricci scalar an additional effective field of the theory 
[24| a choice for these additional variables is given by 



(44) 



where lZ a determines the fluctuations in the Ricci scalar A and 5i a the ones of its momentum A. As for the other 
perturbation variables here we will consider only their scalar part 



A = S 2 V 2 R, $t = S 2 V 2 R 



(45) 



The results of [25[ will guarantee that these new quantities are indeed gauge invariant. The set of variables A, Z, C, A, 3? 
completely characterizes the evolution of the density perturbations in /(A)-gravity once the background is given. 

Using equations (|43|) one can derive the evolution and constraint equations for the variables (|42|) and (|45|) . These 
equations constitute a system of first order partial differential equations [HI, [2l| : 



A m = iu6A m - (1 + w)Z , 



n 

Si 



C 



Rf" 26 
~J> 3~ 



z 



(w - l)(3u? + 2) n 2w0 2 + 3w((j, R + Sp") 1 



2(w + l) f< 



2 2/'/' /'/' I/'. 



RQ 



r 



6(to + l) 



K 



A, 



gr 



5R+ 



-V 2 A m - ^V 2 ^ , 



= 3?- 



1 



-i? A r 



Q + 2R^ r ^J Vt-RZ- 



(3w — 1) /i 



-i? 2 



1/1 

3/" ' /' 
"l8/"7e 18A, 



5 2 e/' s" 2 e 



. f(3) 

A" 



A- 



3 /" 3(w+l) 
_ R 
3 



A 



5 2 e 



./•" 

C + A r 



■R + v 2 n, 

2{w- 1)6 



1 



A, 



6^ 

e 



e/" 



6f 



GROf'f^ - f" (3/ - 2 (9 2 - 3^ fl ) /' + 6A6J" 



4w5 2 6 
3(10 + 1)' 



2 fii 



A, 



25 2 / 
/' 



■51? 



6(/') 2 
25 2 f 0/" - 3A/( 3 ) 



n 



3/' 



A 



(46a) 

(46b) 
(46c) 

(46d) 



(46e) 



together with the constraint 




. fO) f" 
2RQ J — - I 



A 



26/" 2/ 



/' 



^--^V^A^O. (47) 



where we have assumed /"(A) 7^ i.e. we are excluding the GR case. Note that this system is made up of four 
equations, which means that the evolution of every single perturbation variable is determined by a fourth order 
differential equation. This differs from the case of GR, in which these equations are of second order, and it means that 
in general the solutions will have four modes. Hence the perturbations in these theories have a much richer structure 
with respect to the GR ones. In addition to that, the coefficients depend on derivatives of the scale factor of order up 
to four, so that the behavior of the perturbation is much more sensitive to the feature of the background than GR. 

In order to reduce the system above to a set of ordinary differential equations, one defines the eigenfunctions of the 
spatial Laplacc-Beltrami operator: 



V 2 Q = --^Q, 



(48) 
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where k = 2ttS/X is the wave number and Q — 0, and expands every first order quantity in the above equations: 

X(i,x) = ]Tx( fc )(*)Q( fe >(x), (49) 

where ^ stands for both summation over discrete or integration over continuous indices. In this way, one obtains the 
equations describing the k th mode for scalar perturbations in f(R) gravity. They are [l8l[2l|: 



Rf" 



r 



(w- l)(3w + 2)n m 2w& 2 + 3w{^ R + 3p R ) 



2 /' S 2 2 /' /' /' /' + 



2(w + l) f 

r 
r 



\ 2 . f (3) 

i86 V 



6(k; + 1) 



w + 1 



. f(3) 



Sjj(fc) _ ^ Z (fc) 



(3tu - 1) 



P (' 1/ ' + /(4) ^ + B/? /(3) ■ 



lBf'TZ 18A, ; 



s 2 Qf s 2 e 



7? 
A' 



3 /" 3(tu + l) 

/" 3" 

3 



R 



AL fc) - 



5 2 e 



C + A 



2(k; - 1)6 6/i fl 

w + 1 e~ 



6f // 6i?e/'/( 3 ) - 6fc 2 /'7' + /" 3/ - 2 (e 2 - 3m r ) f + 6i?ef 

5? H 

/' 



S 2 



4 ^ 2@ -A« + ^^) 



9(/') 2 
25 2 f Qf" - 3Rf& 



K 



3(w + l) 



3/' 



,£(*). 



e/" 



v 2 ^ 



(50a) 

(50b) 
(50c) 

(50d) 



(50e) 



,/ (3) /* 



= 



/' (/o 2 v J ; /' 5 2 _ 

This system takes a more manageable form if we reduce it to a pair of second order equations: 

AW + AA$ + BA$ = CK {k) + Vll^ , 

f"ii {k) + sn {k} + rii {k) = gA$ + ha$ . 



(51) 



(52a) 
(52b) 
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where 



A 
B 



2 \ n Rf " 
4 „■,:«„»■„», ^ 



C = -(«; + !) 



i-2 111 , . /■« _ W3) 

-V7" 1 + ( /_ 2//n + 2i?e/ ") 7^ _ 2i?e 7r 



(w + i)er 



f = ( e/" + 2i?/< 3 ) 

2 

S*2 J ' ~S 2 



T 



fi f 

3 + 6^ " + ^ V ~ ^ /(3) ~ 9/(3)i? ~ /(4) ^ 



H 



I(3u,-iy» + -^- 

3 1 + W 

(w - l)Rf" 
w + 1 



7 • 



f^R 2 + (p R + n R )f + ^RQf" + Rf" 



(53a) 
(53b) 

(53c) 

(53d) 
(53e) 

(53f) 
(53g) 
(53h) 



In the f(R) = R case these equations reduce to the standard equations for the evolution of the scalar perturbations 
in GR: 



AW 

K {k) = {3w - 1) fx m A^ . 



w 



EL 

s 2 



A {k) = 



(54a) 
(54b) 



If one compares the system (|52p with the equations for the evolution of scalar perturbations for two interacting fluids in 
GR one notices that they have the same structure, i.e., one finds friction terms and source terms due to the interaction 
and the gravitation of the two effective fluids. It is then natural to ask ourselves if this analogy can be useful to better 
understand the physics of these models. The answer is affirmative, but with some very important caveats. First of 

all a more correct way to draw this analogy would be to write the system of equations for A m and A# = 5 and 
analyze their structure rather than using the ones above. Also, as stated in [l8|], one has to be careful in remembering 
that we are dealing with effective fluids and, as such, they might violate some basic constraints that standard fluids 
usually follow (such as the energy conditions) or present subtleties in the definition of their comoving frame. 

However, in spite of these differences one can still use the coefficients of the (A m , Ar) equations to obtain information 
about the interaction between standard matter and the curvature fluid. In fact, the coefficients of these equations 
are found to behave as a ratio of polynomials in the wavenumber and have a non trivial behavior in t. This kind 
of behavior is very different to what is found in a GR-two fluid system. In a photon-baryon system, for example, 
the dissipation terms grow as k 2 and behave like 1/t in time. Therefore, unlike the Thompson scattering in the 
baryon-photon system, the effect of the interaction between matter and non-linear gravitation can influence large and 
small scales alike, depending on the structure of the action. Hence, from the distribution of the structures in the 
observed sky one can deduce constraints on the nature of the theory of gravity. Another important difference with 
GR is that (f5"2"j) are scale dependent for any value of the barotropic factor. This means that whatever the equation 
of state of standard matter, the perturbation solutions will always depend on the scale at which they are calculated. 
This is true even in the special case of dust which in GR is associated to a scale invariant perturbation. 

An additional interesting feature is the meaning of the long and short wavelength limit in (|52[) . As it is clear from 

(|54[) , in GR these limits can be defined by comparing the values of the quantity S 2qi and the ratio between matter 
term in the A m coefficient and O 2 . However, looking at (|52p it is clear that the situation here is more delicate. For 
example, the short wavelengths regime cannot be defined as simply s % &i ^> 1, but S iQ>i has to be bigger than all 
the other quantities appearing in the coefficients B, C and J- divided by 2 . The same reasoning holds for the long 
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wavelengths: g ^ has to be smaller than all the other quantities appearing in the coefficients B, C and T 3 . In this 
sense the structure of (|52|) suggests the presence of a least three different regimes in the evolution of the perturbations. 
Firstly, the "deep super-horizon" regime in which k is effectively zero, an intermediate one (or two depending on the 
value of the barotropic factor w) which is determined by the details of the background and a "deep subhorizon" regime 
in which one has effectively k — > oo. The situation seems to become even more complicated when the fourth order 
gravity action posses dimensional constants (like in the case f(R) = R + aR n ). Since these constants are associated 
with the scales at which the different contributions to the action become dominant, one would expect the introduction 
of additional scales into the theory, i.e., further possible evolution regimes for scalar perturbations. 

V. EXAMPLES 

The considerations drawn above are the most that we can extract from the system of perturbation equations in its 
general form (|52[) . In order to obtain more information we need to apply these equations to some specific models. The 
main issue related to this, is the lack of exact or numerical backgrounds in the most popular models of /(-R)-gravity. 
This is due on a side to the complexity of the equations which makes very difficult their resolution, on the other to the 
difficulty in findings constraints on the values of the parameters which are needed to perform numerical integration. 
The error due to the use of a background which is not a solution of the cosmological equations is very difficult to 
estimate and it could well depend on the detailed of the model considered. For this reason we will investigate here 
some very simple models for which some backgrounds arc known in order to gather as much information as possible 
and test our general conclusions. 

A. The case of i? n -gravity. 

Let us consider the case f{R) = X-R", called also sometimes i?™-gravity, whose action reads 

L = V^g[xR n + 2C m ] , (55) 

where C m is the matter Lagrangian. The above theory constitutes the simplest possible example of fourth-order 
gravity. In this toy model the cosmological equations associated with a Friedmann-Lemaitre-Robertson -Walker 
(FRLW) metric are particularly easy to deal with. The FLRW dynamics of this model has been investigated via a 
complete phase space analysis in Q. This analysis shows that for specific intervals of the parameter n there is a set 
of initial conditions with non-zero measure for which the cosmic histories include a transient decelerated phase which 
evolves towards an accelerated expansion one. This first phase, characterized by 

S = SotW^l , K = , //" = y™r 2n , (56) 
was argued to be suitable for the structure formation to take place. For this reason we will use (|56|) as background. 

1. The long wavelength limit 

Let us start analyzing the evolution of the long wavelength perturbations during the phase (|56|) . In the long 
wavelength limit the wavenumber k is considered to be so small that the wavelength A = I^Sjk associated with it 
is much larger than the Hubble radius. Equation (|48| then implies that all the Laplacians can be neglected and the 
spatial dependence of the perturbation variables can be factored out. In GR, one can prove [26[ that in this limit 
and in spatially flat {K = 0) backgrounds the equations for the variable C reduces to C = i.e. the variable C is 
conserved. From the structure of equations (|46e[) above one can see that in this limit and in spatially flat (K = 0) 
backgrounds the (|50c|l reduces to C — i.e. the variable C is conserved also in the f(R) case. This means that the 
number of equations in (|50p can be reduced to three. 



However, there is a difference between the two limits because in the long wavelength limit one can always set giQ'i ~ 0' while in the 
short wavelength the definition is completely dependent on the values of the coefficients and. as consequences, on the features of the 
background. 
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Substituting the background (1561) in the equations thus obtained and decupling the density perturbation equation 
one obtains: 

(n - l)A m - (n - 1) ( i^- - 5^ ^ + 2? x (n, «>) ^ + P a (n, "0^ + «;) C r - 1 = , (57) 
\ w + 1 J t t z t 

where, Cq is the conserved value for the quantity C and 

2(-9(2(n-l)n + l)w 2 + 6n(n(4n-7) + l)w + 18w + ?i(4n(8n-19) + 33) + 9) , 

Vl{n) = 9(^Ti7 ' (58) 

_ . , 2((2n - l)w - l)(4n - 3(w + l))(3(io + 1) + n(-9io + n(6w + 8) - 13)) 

^(») = 9^rrr ' (59) 



V 3 (n) = — " _ u \ w ~<~ ' _ v / f (60) 



n(21to - 6n(w + 2) + 31) - 18(w + 1) 
This equation admits the general solution 



A m = K^' 1 + K 2 t a + + K 3 t a - - K A t 2 ~W^ , (61) 



where 



1 nw , J{n - 1) (4(3u> + 8) 2 n 3 - 4(3w(18u; + 55) + 152)n 2 + 3(w + l)(87w + 139)n - 81(w + 1)) 
a± = — — H — ± 



2 



2 w + 1 6(n- l)(w + 1) 

(62) 

= 9(w + l) 3 (18(m + 1) + n(-21w + 6n(w + 2) - 31)) 

4 ~ 8{n(6w + 4) - 9(w + 1)) (6(iu + 2)n 3 - (9w + 19)n 2 - 3(w + l)(3w + l)n + 9 (to + l) 2 ) ' [ ' 

Let us now focus on the case of dust (w = 0) . The above solution becomes 

A m = K.t- 1 + K 2 t a +^=» + K 3 t a -^=° - K^t 1 -^ , (64) 



: Q2 
°0 



where 



1 , y/(n - l)(n(32n(8n - 19) + 417) - 81) 

u±\ w =o = --± -r , 65) 

2 6(n — 1) 

K I - 9(rc(12n-31) + 18) 

4|u,=0 ~ 8(4n - 9) (12n 3 - 19n 2 - 3n + 9) " 1 ' 

A graphical representation of the behavior of the exponent of the modes of (|6~4"|) as n changes is given in Figure [TJ This 
solution has many interesting features. For 0.33 < n < 0.71 and 1 < n < 1.32 4 the modes t a± ^=° become oscillatory. 
However, since the real part of the exponents a±\ w =Q is always negative, the oscillation are damped and bound to 
become subdominant at late times. The appearance of this kind of modes is not associated with any peculiar behavior 
of the thermodynamic quantities in the background i.e. none of the energy condition are violated for the values of n 
which are associated with the oscillations. The nature of these oscillations is then an higher order phenomenon. Also, 
for most of the values of n the perturbations grow faster in i? ra -gravity than in GR. In fact only for 1.32 < n < 1.43 
all the modes grow with a rate slower than t 2 ^ 3 . 

Probably the most striking feature of the solutions and (|6ip is that the long wavelength perturbations can grow 
if the Universe is in a state of accelerated expansion (see Figured]). The consequence of this feature is quite impressive 
because it implies, against our (GR-based) intuition, that in R n gravity large scale structures can in principle also be 
formed in accelerating backgrounds. The suppression of perturbations due to the presence of classical forms of Dark 
Energy (DE) is one of the most important sources of constraints on the nature of DE itself. Our example shows that 



4 The values of n presented here and in the following are the result of the resolution of algebraic equations of order greater than two. These 
values have been calculated numerically. This means that the values we will give for the interval n will be necessarily an approximation. 
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if one considers DE as a manifestation of the non-Einstcinian nature of the gravitational interaction on large scales, 
there is the possibility to have an accelerated expanding background that is compatible with the growth of structures. 

In the limit n — > 1 two of the modes of (|61|) reproduce the two classical modes t 2 ^ 3 and i _1 typical of GR, but the 
other two diverge. At first glance this might be surprising, but it does not represent a real pathology of the model. 
In fact, when n = 1, equation (|57[) reduces to a first order differential equation with a forcing term, which presents 
a two modes solution. This implies that in this case the additional modes in the solution (|61j) can be discarded and 
GR is recovered. 

Solving for the other variables we can also obtain the solution for the other scalars: 

si 

K = K 5 t%Ti- 3 + K 6t P+ + Krt 13 - - Kg -^t~ ( 67 ) 
ft = Kgt^- 1 + K w t~<+ + K X1 P- ~ K^t-^TT)- 1 (68) 

where 



P± = a± - 2 , (69) 
7± - a±-3, (70) 

and the constants K§, ..Kyi are all functions of K±, ..K4. These expressions are rather complicated and will not be 
given here. It is interesting that these quantities have an oscillatory behavior for the same values of n for which A m 
is oscillating. Also for these quantities the oscillating modes are always decreasing. 




2. Short wavelengths. 



The analysis of the properties of the small scale perturbations is more complicated than long wavelength ones. This 
because, if, as often happens, it is not possible to obtain exact solutions in this regime one is bound to find other 
means the relevant physical information. 

A way of doing that, based on a statistical treatment of the behavior of the perturbation at different scales, is 
to define some characteristic quantities which help extract physical content from these solutions. One of the most 
important of these quantities is the power spectrum Px{k) of a perturbation variable X. Mathematically Px{k) is 
defined as the variance of the amplitudes of the Fourier transform of A at a given value of k [27j : 



(A(kx)A(k 2 )) = P x (k 1 )6(k 1 +k 2 ) , 



(71) 
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where arc two wavevectors characterizing two Fourier components of X and P(ki) = P(ki) because of isotropy 
in the distribution of the perturbations. The power spectrum Px(k) tells us how the fluctuations associated to X 
depend on the wavenumber at a specific time and carries information about the amplitude of these perturbations (but 
not on their spatial structure). Its features are directly connected with the structure of the equations. For example, 
(|54[) shows that in the case of GR and pure dust the matter fluctuations equation does not depend on k , so the matter 
power spectrum can be considered constant (i.e. scale invariant) [28| . 

A spectrum of one of the perturbation variables is a powerful tool for comparing the predictions of the system 
(|52p with observations and is able to reveal a great deal of information on the physics of scalar perturbations. In 
what follows we will focus on the matter perturbation spectrum Pa(&) only, because this spectrum is the one that 
can be directly measured via large scale surveys. Our spectra will be normalized in such a way that it is unity on 
super-horizon scales. This normalization can then be scaled with current observations of the power spectrum on large 
scales (see [29] for the latest constraints). We will also analyze the time variation of this spectrum in order to obtain 
information on the ways in which the perturbations evolve in time on different scales. 

Let us then consider Pa(&) for (f52|) at t = S/Sq = 1 (Figures 2(a) 2(b) [3]). As mentioned above the fc-structure 
of Equations (j52|) suggests that in fourth order gravity there exist at least three different growth regimes of the 
perturbations. This is confirmed by our results for this example. In particular, in the case of dust and for any values 
of the remaining parameters we have that: (i) on very large scales the spectrum goes like GR i.e. it is scale invariant; 
(ii) as k becomes bigger the scale invariance is broken and oscillations in the spectrum appear; (iii) for even larger k 
the spectrum becomes again scale invariant. However, on these scales the spectrum can contain either an excess or 
deficit of power depending on the value of n. In particular for n ~ 1 + small scales have more power than large scales, 
but, as one moves towards larger values of n, the small scale modes are suppressed. For < n < 1, instead the drop 
in power seems to decrease as one moves from n = towards 1~ and we see a sudden increase for n ~ 1 _ . It is worth 
noting the case n « 0.8 for which there is basically no difference in power between large and small scales and there 
are no significant oscillations in the spectrum. 

Further indication of the link between the k structure of (|52|) and the different regimes of the matter power spectrum 
can be seen if one analyzes the power spectrum in a radiation dominated era. In this case the perturbations equations 
contain an additional fc-term which is not present in the dust case. This means that one would expect four different 
regimes, rather than three. In Figure [4] we have plotted the matter power spectrum of i?"-gravity in the case of 
radiation and n — 10 (this value of n is chosen only for convenience), and as expected one can recognize four different 
regimes. Note that we obtain the same k scaling as in GR on small scales. 

Finally, further information on the dynamics of the matter perturbations can be obtained examining the time 
evolution of the power spectrum. In Figure [5] we give the power spectrum for n = 1.4 at different times. One can see 
that, as the universe expands, the small scale part of the spectrum is more and more suppressed and oscillations start 
to form. On the other hand the large scales do not seem to be evolving, which might appear in contrast with what 
mentioned above. However, this is a byproduct of the normalization: for clarity we have normalized the spectrum in 
such a way that every curve has the same power in long wavelength limit. 

The features of the spectra that we have derived can be then interpreted in terms of the interaction between the 
curvature fluid and standard matter plj . On very large and very small scales, the coefficients of equations (|52[) become 
independent from k so that the evolution of the perturbations does not change with the scale and the power spectrum 
is scale invariant. On intermediate scales the interaction between the two fluids is maximized and the curvature fluid 
acts as a relativistic component whose pressure is responsible for the oscillations and the dissipation of the small scale 
perturbations in the same way in which the photons operate in a baryon-photon system 5 . The result is a considerable 
loss of power for a relatively small variation of the parameter n. For example, in the case n = 1.4 the difference in 
power between the two scale invariant parts of the spectrum for n = 1.1 is of one order of magnitude while for n = 1.6 
is about ten orders of magnitude. 

Probably the most important consequence of the form of the spectrum presented above is the fact that the effect 
of these type of fourth order corrections is evident only for a special range of scales, while the rest of the spectrum 
has the same k dependence of GR (but different amplitude) . This implies that we have a spectrum that both satisfies 
the requirement for scale invariance and has distinct features that one could in principle detect, by combining future 
Cosmic Microwave Background (CMB) and large scale surveys (LSS) (30l. [3l|. 



5 This suggests the following interesting interpretation for the perturbation variables 1Z and 5R. These quantities can be thought to 
represent the modes associated with the contribution of the additional scalar degree of freedom typical of f(R)- gravity. In this sense 
the spectrum can be explained physically as a consequence of the interaction between these scalar modes, which function as a scalar 
gravitational wave, and standard matter. 
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(a)Plot of Pa(&) at t = 1 for /{"-gravity and n > 1. Note that the spectrum is composed of three parts 
corresponding to three different evolution regimes for the perturbations. 
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(b)Plot of the P&(k) as a function of k for i? n -gravity at t = 1 and 
< n < 1. The spectra for n = 0.6 and n = 0.5 only approach the scale 
invariant plateau at extremely high k when compared to the other curves. 

Note the behavior of the spectrum for n m 0.8, differently from all the 
other cases, there is basically no loss of power in the spectrum at large k. 

FIG. 2: Matter power spectra for Ji n -gravity (from [2l[) 
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FIG. 3: Detail of the plot of Pa(&) for i? n -gravity at r — 1. As expected, for these values of n we find the presence of the three 
regimes mentioned in the text. Note also that for n ~ 1 + the small scales are characterized by an excess of power [2l]]. Such 
features are compatible with some of the results found in [{J. 
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FIG. 4: Pa(&) for i? n -gravity at r = 1 in the case of radiation for GR and i?™-gravity with n = 10. As expected in this last 
plot we find four different regimes instead of the three of the dust case: a first regime for k — > which is scale invariant; a 
second and a third regime which correspond to the two different slopes between k w 10 -2 and k ~ 10 1//2 and a fourth regime 
which has the same slope of the GR plot [2l| . 
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FIG. 5: Evolution of -Pa(&) for i?"-gravity for n = 1.4. The spectrum has been normalized in such a way that the curves 
coincide at large scales. Note how, as time passes, small scale perturbations are dissipated and oscillation appear [2lj . 



B. The case f(R) = R + aR n 

Let us now consider an f(R) theory characterized by the following action, 

L = V^fJ [R + aR n + C M ] ■ (72) 

This th eory has gained much popularity as a fourth order gravity model within the context of both inflation and dark 
energy p33M36j. 

Unlike i?"-gravity, (|72p includes explicitly the Hilbert-Einstein term, so that one can consider it as the result of an 
additive correction to GR. As a consequence this model posses an additional physical scale which is related, in our 
units, to the coupling constant a by the relation L = a 2 ^™ _1 \ making (fT2"j) the simplest fourth order gravity theory 
with an additional scale for the gravitational interaction. 

In the following we will take a, which in our units is the ratio between the coupling constant of the fourth order 
corrections, to be positive definite. Of course, we expect this model to behave much in the same way as the model 
discussed in the previous section for a — > oo, and to recover GR for a = 0. This also means that any new feature in 
this model will emerge for intermediate values of the coupling. 

The f(R) = R + aR n model has been analyzed at the level of the background using many different approaches 
(see for example [32T[36l |). but probably some of the most interesting results for cosmology have been found using 
the dynamical system approach H, The dynamical systems analysis proved that this class of models has, like 
i?™-gravity, an unstable fixed point associated with the Friedmann-like solution 6 a = t 2n /H 1 + w ) _ 

One can then analyze the behavior of the perturbations of this background in f(R) = R + aR n and compare it 
with the one found in i?"-gravity. However, dealing with this model is considerably more complicated than dealing 



Unlike the previous example, the structure of the phase space for this model is not well known. This means that, although in this theory 
one has a fixed point that resembles the one of i? n -gravity, it is not obvious that it plays the same role. In addition, this background 
is not, in general, a physical solution of the cosmological equations Q. This means that there are cosmic histories in which the general 
integral of these equations approximates the behavior we consider, but it will never be exactly the same. We choose to treat this as a 
further approximation in our investigation. 
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with i? n -gravity. For example, one is unable to find exact solutions for the perturbation equations, even in the long 
wavelength limit. For this reason in the following we will focus directly on the features of the power spectrum which 
are shown for various values of the parameters a and n in Figures El As one can clearly see, these plots resemble 
the ones we have derived in the previous example. There are, of course, differences in the position of the oscillations 
and the amount of the power drop, but one finds again three different regimes in the case of dust and two of them 
(k — > and k — > oo ) correspond to scale invariance. Particularly interesting is the fact that in principle the values 
of a and n can be fine tuned in such a way to obtain a spectrum in which the small scales have the same power as 
the large ones. In a situation like this most of the spectrum would be scale invariant and all the deviations would be 
concentrated around a specific scale. 

The time evolution of these spectra also reveals some interesting insights into the dynamics of the matter fluctuations. 
As usual for large values of a the evolution is very similar to the one obtained for i?™-gravity as it is shown in Figure 



7(a) However, when the value of the coupling changes the behavior of the perturbations can change dramatically. An 
example is given in figure [7(b) | representing the power spectrum of the model (n = 1.4, a = 0.01) in which the small 
scale perturbations are first dissipated and successively start to grow again. This means that in principle one could 
choose n and a, such that for example the small scale perturbation grow at different rates at different times. This 
property could be useful in the resolution of open problems in GR structure formation, like the cosmological dark 
matter or the excess of dwarf red galaxies. 

Therefore, in spite of the presence of an additional scale, the power spectrum in this class of model seems to 
preserve most of the main structure of the one in i?™-gravity. This implies that all the considerations made in the 
previous section concerning the physical mechanisms behind the form of the spectrum can be made also in this case. 
At this point the natural question that arises is about the generality of this result. Is it possible that we have found 
a characteristic feature of /(i?)-gravity? Unfortunately, due to the lack of exact background solutions (many of the 
more complicated models of f(R) do not admit power law solutions) and the difficulties related with the numerical 
resolution of the perturbation equations, one is unable to answer directly to this question. However, the important 
point is that if, as we have argued in the previous sections, the structure of the power spectrum is determined only by 
the /c-structure of the perturbation equations, one can make the conjecture that the features we have observed in our 
two simple models are indeed common to all the f(R) models. If this is confirmed then we would have found a crucial 
method to test the existence of higher order corrections and it would be a powerful constrain for the parameters of 
any f(R) model. Another issue concerns the relevance of the form of the background in emergence of the features 
listed above. Can specific choices of background influence the presence of this characteristic feature? In order to 
answer this question we will consider another simple example: a de Sitter Universe. 



C. Perturbations of the de Sitter spacetime in a general f(R) -gravity. 

Let us consider now the cosmological perturbations around a de Sitter background in /(i?)-gravity. The presence 
of a de Sitter background in /(i?)-gravity is one of the most important features of these theories because it has the 
potential to model both inflation and dark energy [3^, [3j|. In fact, it has been proven that a viable /(i?)-gravity 
model unifying inflation and late time acceleration in the form of double de Sitter solution can be always constructed 
numerically |40j. As we will see, however, such backgrounds are not suitable for structure formation, because matter 
is dissipated very quickly. Notwithstanding this physical issue, the peculiar properties of this metric allow us to go 
deeper in the understanding of the perturbation equations. 

Let us consider a Universe in which the background is given by a de Sitter spacetime characterized by a scale factor 
S = Soe* 9 * and vacuum (fi m = 0). Substituting in the cosmological equations it is easy to show that f3 has to satisfy 
the equation 

18/3 2 /o-/o = 0, (73) 

where /q = f'(Ro), fo = f(Ro) and Ro = 12/3 2 . The (f73"f also implies that not all the f{R) theories of gravity admit 
de Sitters solutions, consistently with what one finds in 

Let us now consider a perturbation of this spacetime in which a fluid, constituted for instance by standard matter, is 
present. According to what has been said in the previous sections, this fluid will be described by first-order quantities. 
We will also assume that the fluid is actually barotropic in its rest frame i.e. its equation of state is p m = w/i m . 
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(a)Plot of Pa(&) for n > 1, a = 10 and various values of k 



(b)Plot of Pa(^) for n > 1, a = 1 and various values of fc 
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(c)Plot of Pa(^) for n > 1, a = 0.1 and various values of k (d)Plot of P/±(k) for n > 1, a = 0.01 and various values of fc 



FIG. 6: Plot of Pa (it) as a function of k for P + o?P"-gravity at t = 1 for n > 1 [21 
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(a)Plot of Pa(^) for n = 1.4, a = 10 and evaluated at (b)Plot of P/\(k) for n = 1.4, a = 0.01 and evaluated at 

various values of r various values of t 



FIG. 7: The time evolution of the matter power spectrum in R + Q_R™-gravity for n = 1.4 and differing values of a. Note the 
drastically differing vertical scales in the plots. Additionally, note the rise in power at small scales at late times in the case of 
n = 1.4 and a = 0.01 EOf. 



Choosing a set of observers comoving with it 7 , the perturbations equations turn out to be [41[ 



A = -eo 

K = U, 

ft = -Gft-\ 
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(74a) 
(74b) 

(74c) 



where, as usual, we have assumed f"(R) ^ i.e. we are excluding the GR case. Note that these equations cannot be 
obtained plugging the background we defined directly from (|46p . This because the fact that the matter thermodynam- 
ical variables are of first order changes the structure of linearized 1+3 equations, which, in turn, leads to changes in 
the differential structure of the equations rather than only in the coefficients. Performing the harmonic decomposition 
and substituting (|74b|) in (|74cj) we obtain j4l| 



= -3/3(iu ■ 



l)A( fe ) , 



12/3 2 + ^ 
Jo 



(1 -3w)AW 
3/£ 



0. 



(75a) 
(75b) 



In this system the equation for A is scale invariant and, as expected, matter perturbations are exponentially sup- 
pressed with a time constant which depends on w and the time constant of the de Sitter solution. The Ricci scalar 
perturbations, instead, are governed by a second-order equation which is forced by the matter term. 
In the long wavelength limit k = the above equations yield the general solutions 



A = A e 



-3t/3(l+w) 



Ac 



-3(™+l) 



(76) 
(77) 



7 Since the definition of the fluid flow u a is made at the level of the perturbed Universe this choice is legitimate. In addition to that, the 
dc Sitter solution is frame invariant so any choice of frame in the background is equivalent. 
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where 



a±=-3/3±y25/? 2 -^|. (78) 

Here IZo.i and Ao are constants of integration and we have dropped the apex "(0)" to make the notation lighter. It 
is plain from this solution that, in a de Sitter background, standard matter is clearly made homogeneous, but this is 
not the case for the perturbation of the Ricci curvature. If one considers 1Z as representing the scalar gravitational 
waves normally associated to the scalar degree of freedom of this type of theories, one can see that, depending on the 
form of the function /, this kind of perturbation is able to grow. In addition, if we imagine our /(i?)-model to be an 
inflationary one, we can see that the analysis of the scalar waves would constitute a direct and purely classical test of 
the nature of the gravitational interaction, based on the gravitational wave relics of the inflationary era. 

The form of the exponents of the modes the solutions above for some popular models of /(i?)-gravity [22|, HH, |39| 
are given in Table |U 

It is easy now to derive the properties of the power spectrum for the matter perturbations. It is clear from (|75al) 
that their spectrum it is scale invariant and remains constantly scale invariant at all times. Comparing these results 
with the ones of the previous section one realizes that f(R) de Sitter inflation is very different from an f(R) power 
law inflation. This result is similar to what happens in GR when de Sitter and power law inflation are compared: also 
in this case power law inflation is accompanied by a loss of power in the infrared part of the spectrum j44| . 

The fact that with a de Sitter background the power spectrum of the theory is flat seems to be in contrast with 
the claim of the presence of characteristic signature made in the previous section. However one must remember that 
the reason why this happens is the very special properties of the de Sitter metric and the fact that in our background 
there is no matter. This eliminates all the k dependence from the equation of the matter perturbations giving rise 
to a scale invariant spectrum. Such result reminds us that, independently from the form of /(-R), a specific choice 
of the background can change deeply the system (|50p and conditions the appearance of our characteristic signature. 
Therefore one can conclude that, in the appearance of the characteristic signature in the power spectrum the structure 
of the perturbation equations due to the background is more important than the additional scales contained in the 
function f(R). This result shows again how the choice of a correct background is crucial to make any reliable prediction 
in perturbation theory. 



VI. DISCUSSION AND CONCLUSIONS 



In this paper I have reviewed the construction of a covariant and gauge invariant formalism that allows the descrip- 
tion the perturbations around any cosmological background for a generic f(R) theory of gravity. Using this technique 
I have been able to derive, in a frame comoving with matter, a system of equations that describes the evolution of 
the scalar perturbations which are traditionally associated with structure formation. The covariance of the formalism 
guarantees that we are able to write these equations in any other frame. Differently from the case of GR, this system 
appears to be of order four rather two, and the coefficient of the equations depends on derivatives up to four of the scale 
factor. Also, in their harmonic decomposed form, the equations present a non trivial dependence on the wavenumber 
k, so that perturbations of dust can depend on the scale in spite of the absence of fluid pressure. Such features are 
originated by the additional Laplacians appearing in the equations which, in turn, are a direct consequence of the 
fourth order terms in the gravitational field equations. 

The application of these equations to some simple /(i?)-theories helps to further clarify their properties. We have 
seen that the behavior of the perturbations in these theories can be very different from the one in GR. For example, 
in i? ra -gravity long wavelength perturbations can grow in a power law inflationary scenario. Also we learned how 
important are the features of the background in the determination of the behavior of the perturbations. For example, 
form the analysis of the de Sitter backgrounds one learns that the specific fc-structure of the equations can be changed 
by the use of backgrounds with special properties. Our results also help in the understanding of more general features 
of perturbations in /(i?)-gravity. The comparison of the results in R n and R + aR n on the same background, shows 
us that the presence of an additional universal constant in the action (and the associated lenghtscale) , although 
introducing changes in the dynamics, has little influence on the power spectrum. In this sense the structure of the 
power spectrum, should preserve some basic characteristics, like the presence of at least three regimes of which two 
are scale invariant (but not necessarily at the same power). 

Can we speak then of "characteristic signature" of /(i?)-gravity? And in what sense? The information that we 
have discovered so far seems to show that the introduction of additional scales i.e. additive terms in the function / 
does not affect the basic features of the matter power spectrum, because they are only related to the fc-structure of 
the perturbation equations. However, such structure can be modified by the specific features of the background. In 
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TABLE I: Some of the values of ft and the exponents of the modes of the scalar fluctuation solutions for various popular f(R)- 
gravity models in pure de Sitter backgrounds [4l| . For the more complex forms of f(R) the implicit equations to be solved 
in order to find the parameters have been given. Of special interest are the models f(R) = i+^j;* an d their generalizations, 
which can provide a unique theoretical framework for early time inflation and late time acceleration [38|, [3^ (the first unified 
models of this type were proposed in |43| ) . 
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this sense one can conclude that a class of background which do not change /c-structure of the equation are bound 
to generate for any (non pathological) f(R) theory spectra with similar features in terms of regimes and behavior at 
large and small scales. This similarity is what we call a "characteristic signature" of these models. 

On a more observational point of view the important consequence of this conclusion is that a deviation from the 
scale invariance of the power spectrum can be interpreted as caused by non Einstcnian gravitational interaction. 
Therefore the large scale structure surveys can be used as powerful tests for the gravitational interaction. 

The same hold in the case of inflationary models based on f(R) gravity. As in the case of GR, power law inflation 
is different from de Sitter inflation, and the detection of scalar gravitational waves, could not only be a footprint of 
non Einsteinian gravity at work, but also could allow us to distinguish between power law inflation and de Sitter one, 
and even give information on what form of f(R) gravity is most likely to be at work. 

All these results look very promising. In spite of the many difference between /(i?)-gravity and GR, the observable 
seem not completely in contrast with the observations. Of course the next step is now to try to refine the rough 
models we have used up to now to obtained better results. These will help shedding more light on this fascinating 
topic. 
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